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Ahoccncc 

A |>roct'(liiro Iti pruHaiiLud for tlio Jictlvo con- 
trol of n iipitmlnp, flexible tipncocraft. Rucb n 
iiyiilctii oxlilbllii |iyn))ii' 0 |ilo errooLii. Tlui di'iilpii 
of llio controller In booed on modnl doconipooltlon 
of till* p.yroncoplc ityotcin, Tlilfi modnl dacmipllnij 
procedure loiido to a control muclinninm Implemented 
Jn modular formi wlilcb repreoentn n dlntlnct com- 
puLAtlonnl ndvnntnp.o over the control of tlio 
coupled oyoteni. IloitlBn proccdureo ore deraonatmted 
for two typoii of control nlcorltbraSi linear and 
nonlinear. Tito flrut representa claaalcnl linear 
feedback approacb and tbc second represmts an up-; 
plication of on-off control, both types made fea- 
sible by the modal decomposition sciicmc. 

I, Introduction 

Ao spacecraft structures Increase in size, 
wclglit limitations demand that various substruc- 
I lures be made as light ns possible, which In turn 
I requires that they be highly flexible. On the 
other hand, greater pointing accuracy necessitates 
finer attitude control, which can be achieved only 
lltrough active control of the spacecraft. In con- 
I trolling a flexible spacecraft, the problem of sl- 
' mulnting the flexibility Is often critical, as the 
I number of degrees of freedom of the simulation can 
I become so large as to render various mathematical 
1 techniques unConslblo, Of course, quite often 
; proper modelling of the spacecraft can result In 
J a model possessing relatively few degrees of free- 
! dom and yet retaining all the essential dynamic 
chntnctoristlcs of the system. Even then some 
truncation may be necessary. In the case of con- 
trol of nonrotnting spacecraft, It is common prac- 
tice to use the system natural modes to decouple 
the system and control only a limited number of 
lower modes. This procedure can be implemented 
with relative ease because natural modes of nonro- 
tating structures con be readily computed. Recent 
advances in the analysis of gyroscopic systems, 
however, makes o modal approach possible also for 
rotating spacecraft. Following is a brief litera- 
ture survey of related work. 

In an attempt to control a flexible space boos- 
ter, Covarter (Ref. 1) presents a procedure where- 
by the response can be represented in terms of the 
rigid-body modes and the bending modes of tbc mis- 
sile. This approach permits a dosetiption of the 
sy.Htem In terms of transfer functions for the un- 
' coupled sy.stcm. When the spacecraft Is spinning, 

I or when it possesses spinning pnrt.s, the classical 
! modal dc'comiiosltlon is no longer possible because 
I the modal matrix will not diagonalize the gyro- 
scopic matrix. Kuo et ul (Ref. 2) liavc presented 
a technique for the design of a digital controller 
for spinning flexible spacecraft using a redesign 
of a preliminary contlnuoee-data control system. 


tio attempt wna mndo In Hof. 2 to uncouple tlio nyii- 
tem cquatlmin. Control of fluxlble spacecraft by 
modnl nyntUoHlii la dliicussud' by I’oaiaort (Kaf. 3), 
but the procedure proiiiMitod lii valid only for non- 
gyroMCupIc ayatnian. In fact, the iiinthematlcnl 
model considered, used first In Kef. >i, consists 
nf three disks mounted on s flexible simfc and ro- 
tating about .1 common symmetry nxln. Such a model 
doca not exhibit tho gyroscopic effect typical of 
a spliming flexible structure copnblc of nutation. 

This paper develops o method for the design of 
a controller based on modal dccompoultion of spin- 
ning structures developed In Refs. 5 and 6. For 
blgh-ordcr systems, this approach offers substan- 
, tlal computational advantages. In the first place, 
the modal decoupling procedure leads to a control 
I mochanloffl which eon be Implemented in modular form, 
i Moreover, one can use decoupled dynamics to design 
an observer also. Following decoupling, each con- 
trol-group is governed by a set of two first-order 
differential equations with a skew symmetric ma- 
trix of coefficients. These sets of equations can 
be integrated readily, thus permitting Independent 
control of spacecraft modes. Design procedures 
are demonstrated for two types of control olgo- 
ricluns, linear and nonlinear. The first repre- 
sents a classical linear feedback approach In the 
form of proportional control and the second repre- 
sents an application of on-off control, both types 
mado feasible by tbc modal dccompoaltlon scheme. 

I 

2. Klncmotlcal Considerations 

Let us consider a general spacecraft consisting 
of a central body with an arbitrary number of ap- 
pendages. The central body will be referred to as 
the "platform" and it can be rigid or clastic. 
Quantities pertaining to the platform will be des- 
ignated by the subscript P. Tlio appendages can be 
of three types: rigid and rotating relative to 

the platform, clastic and nonrotatlng relative to 
the platform, and clastic and rotating relative to 
tbc platform. Quantities pertaining to the types 
of appendages listed will be denoted by the sub- 
i scripts R, E, and A, respectively. An example of 

! the first type is a rigid rotor, examples of the 

second are flexible solar panels or flexible an- 
tennas cantilevered from the platform, and an 11- 
lustrntlon of the third is a flexible rotor. 

Clearly, there can be more thnn one appendage of 
a given type. We shall confine our discussion to 
one of c.nch type, however, with a summntlon Implied 
over appendages of the same type. 

To dcscrlbo the motion of the spacecraft. It 
will prove convenient to Introduce various sets of 
axes. In the first place, we wish to identify an 
Inertial system of axes XYZ with the origin at a 
point 0, Then, wo shall identify a system of axes 
xpypZp with tlic origin at the center of mass P of 
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iilii> t>l»>rnnii nml I'o Inc Id lint will) the pi'JiirL|>itl 
iiixcii III III!' imili'Iiiniicil iilatliu'iii. liinlliirJyi wo j 
(nlinJl coiiiililoi' II net of nxoo nttiiclicil to 

;ilm cnitir nml with the orl);lii lit thu miiiiii cantor K 
Hif till' rotor. The motloii of nn ciiintJc membor iion- 
liotiiiliijt I'l-IiiLivu to the iiliitforin aim bo dciicrlbcd 
(by nicnnii of n oot K|iy|f/.|» nttilcliod to the moiabor In 
iiiiidcl'nnncd iiuiLo niiil' with the origin II nt ihn point 
lof iittiichniouL of the moiiibor. In tlio Hiimu mannuri 
•we cmi dorinc n iiiit of iixun Xy\y/^Zj\ with the origin 
int A mill L'liliicldlin; with ii (tlveii'iiut uf iixiiii 111 i 
itho riUntlng oliintlc mciiibcr when In uiidoforniud 
intate, Tim n|mcccrnft mid the various seta of nxes 
litre iiliuwn ill Klg. 1. 

The pnnltlon of the point P relatlvo to tlio in- | 
crtJnl Gpnee Is given hy the riidlna vector Kop 
ifrom 0 to 1'. The rotation of iixoa Xpyp;;|> relative 
jin nxen XYZ Ja (tlvcn by tlic angular velocity voc- 
jtur (lj>. In addition, any point of thu platform ■ < 
lean undergo clnatlc motion relatlvo to XpypZp. I 
iTor iilmpllclcyt however, we ahall nosumc time tlio 
(platform la rigid. The poaitions of the Intercon- ' 
inerting pointa R, E, and A relntivo to P are de- 
moted by the radiua vcctora Rpp, Rpg, and Rp^i res- 
jpoctivcly, nnd the rotatlons'of axea xpy((Zg and 
|j:,\yA?.A relative to xpypzp are denoted by tlio Piigti- 
jlinr velocity vectors ui[( and respeeri ely. It 
jfoliowa that the poaltlona of R, E, jc.d A relative 
ito 0 are H^p + Rp|j, R{jp + Rpi?, nnd R^p + Rp^, res- 
pectively. The position of an arbitrary point in 
the platform relative to P la given by the radius 
vector rp, Similnrly, the poaltion of an arbitrary 
(l>oint in a spinning rotor relative to R la given 
(by rg. On tlio other hand, points in the flexible 
lappcndagca E and A are described by Jg + ug and 
r. -»• hai where rg and ja are nominal positiono of 
Itlie points when the appendages ore undeformed nnd 
lug mid ua are elastic displacements. It follows 
that the absolute positions of arbitrary points In 


the various 

spsceernft members arc 
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lit should be pointed out that Rqp Is generally 
(given in terras of inertial components, rp, Rpg, 
iRpg, and RpA In Carma of components along axes 
‘Xpypzp, Cg in terms of components along XgygZg, 
[rg T iig in terras of components along xgygzg, and 
Ha yA terms of components along x^yAZA- 

The angular velocity vectors of axes XgygZg, 
XgygZg, nnd Xy^yA^A relative to the inertial space 


ere 

-K ” -P -R ’ ?E " 5p * -A ** “P 

irespcctlvely . Note that Up In teems of components 
laJong I’vy, and iii[{ nnd ray; are in terms of components 
lalong Xpypzp and x^y^^z,;, respectively. Tills per- 
imt tn us to calculate absolute velocity vectors fur 
arbitrary points in the various ceitbcrs in the form 


yp 


YoP 


(3a) 




Hop + Op " Opp + 


Sr 


(3b) 


" Sop Op ** (Bp)-; tn yg) y|i * 

. t 

• V. - v„„ + n„ X U + n X <r. + u.) + u. (3d) 

I -A “OP “P “PA “A “A -A “A 

whore \jg nnd arc velocities of the points In 
<(noiitlon rolntlvQ to the moving fmincs. Wu note 
' once ngnln Hint tlio various terms In R,|s. (3) are 
' In terms of diffsrent (lotii ,of nxos, | 

1 

It will prove Convenient to work wltli velocity 
components In terms of member nxes. Pnr example, 
wo shall express Op nnd vp In terms of componenfo 
along axes xpypZp, ate. To this and It In more 
nntiirnl tn work with matrix notation, wlilcli noeuii- 
sltntcs the Introduction of Llic matrix form uf the 
vector cross product. Hence, let uo define the 
following skew symmetric matrices 
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0 -u 


L“y 


(4) 


Then the cross products SJxr,nxujrx(l, ssid 
u X Q have the matrix counterparts Hr, fit, rH, 
nnd Qn, respectively. In addition, we must intro- 
duce the matrices of direction cosines between var- 
ious systems of axes. For example, the matrix of 
direction cosines between axes XYE and axes xpypZp 
will be denoted by Lqp, so that 


(Xp yp Xp] 


Lop[X Y Z) 


(5) 


' Similarly, the matrix of direction cosines between 
j nxos XpypZp nnd axes Xgygzg will be denoted by Lpg, 
'etc. With this notation, Eejs. (2) can be replaced 
|by 


I ?R “ ^PR-P ’ ?E “ ’'PE^P' 5a " ^PA^'p + «A 

I and Eqs. (3) can be rcplar.od by 
' v„ 


■^R 


^OP-OP “ ’^P?P 

^PR^OP'OP ~ ^PR^PR^p " '^rOr 


.(7) 


,^;e “ SeSpV - Se’^fe^p - "e^’'pe?p + Be 

!ya - h’A^opYop - ^pa^a5p - ^^A “a^5a “a 


3 . Kinetic Energ y, Potential Enoriiv and Hon- 
conservatlvo Virtual Work 

The kinetic energy of any member can be written 
in the gencr,'l form 


J m 




if 

* m_ 




+ I j BeV"e + 1 j yIya‘J"’a 


( 8 ) 


JiriB. (7i, iliu kliioUc onocKX bccomtii "j 

T -|»n vjj, 4| 5? J Sj, + I yj y,^ + | yj 

I 

^ \ % -•■ nj I.;^ Ja *^A + I L 3 **“E 

Big 

■* “l j ^A ‘*"'a " ^OP *'OP^'“e Se ’*' ®A ^A 
”:a 

M T * 

■*■ ^^’JE ^EG ^PE ■*' ^PA *'aG ^PA^5p 

"■ -OP ^OP ^PA *^AC -A Yop ^OP^^'PE Pe ^ ^PA 5a^ 


•+ 5^ "^EG Se ■*■ *Va ^PA "^AG ^PA^5p 

i ^A 'AG ^A - 05e(?e^?eEe ^ ^A^-JaEa^ 
V -yJ'iA 


^T T 
*1 * ■ 


■whETO 311 1 b tlio total mass of the apncccrnlt, VQp 1 b 
tilia unaBnitude of vop, and J Is tUo total Incttla tna- 
itJj: of tiio entire spncocraft In deformed state n- 
})DU»; P ±31 tcimis of components alone xpypzp. More- 

CVCi 


”r 


<10a) 


I 


‘»I+^X'a* V“"a 


EG 


AG 


j (ig + Cg)dmg 

J ^'a “A^‘’^ 


<10b) 


nniioclntcd wJtli tlie olniitlc appenilannii and Dj? and 
I), arc tlio domains of ^uxtimnion of thcou appciidaueo. 
Tno dcnsltlou and depend on spatial durlva- 
tivuB of tlio componentu of the elastic dlsplacc- 
niont vectors n^ I’A* mi'pectlvely. Wo.oliall 
not p,lve explicit cxproualonti for V(.; and at 
thin polntt but return to this subject In the nest 
section. I 


I Tlic potcntlsl energy esn be used to derive tlie 
conservative forcen nctlnp on Uie system, Ln nd- 
dltlon, there can bo noncontiorvatlve forces pres- 
ent, Such forces esn arise from various sources 
Bucli as solar radiation pressure, meteorite im- 
pact, etc. LcttiiiB f be the nonconuarvativc force 
vector per unit area at a given point on the sur- 
face S of the spacecraft and 6K the virtual dls- 
plnccmcnt of that point, the nonconsorvatlvo vir- 
tual work for tlio entire spacecraft can bo wrlt- 
I ton In the form 



fp-‘V=p + Is £»"V=« 

R 

Ce'^V'c^Is £a-‘»»«a 

A 


( 12 ) 


wliccc £Ep, 6Rp, 5Rg, and £Ra can be obtained from 
Eqs, (ij. Note that concentrated forces can be 
treated os distributed by using spatial delta func- 
tions. 


The kinetic energy, potential energy, and vir- 
tual work can be used In conjunction witli Lagrange's 
equations to derive the system equations of motion. 
This system of equations is of the hybrid type, 
l.c. , some of tlie equations arc ordinary differen- 
tial equations and tlio balance are partial differ- 
ential equations. Tlie flr3t arc associated with 
the rigid body motions of the spacecraft whereas 
the second arc associated with the clastic dis- 
placements. It will prove most convenient, how- 
ever, to work with a completely discrete system, 
which requires cite transformation of the partial 
differential equations into sets of ordinary dif- 
ferential equations. This will be done in the 
next section. 


il£ " J ^ ‘‘"’e * Ua " 1 «A '’“a 

JiOg 

i‘z " J ?^^E ‘'"’e 

“e 

(lod) 

i?A " j ^>^A ‘“"a 

’'^A 

Gcncrnlly, the potential energy la of two types, 
namely, gravitational and clastic. Because of the | 
high altitude of gcnsynchi'onous satellites, the i 
dlflerontlal gravity effect Is negligibly small, 
so that the potential energy will he assumed to be ' 
■cntixoly duo to flexibility. Ho shall express the , 
•potential energy in the form 


V » 



dD^-1 



do. 

A 


ai) 


whtsrc Vj, and Vj^ are potential enorgy_dcnBitlc^ 


^ . System Discretization 

To cllminntc the spatial dependence from the 
formulation, let us assume thot the displacement 
vectors ug and can be written in tlio form 


^JE 


E^E 




A^A 


(13) 


where 4>g mid arc rectangular matrices of space- 
dependent admissible functions and Cg and arc 
tlmc-dcpcndcnt vectors of goncroll'zcd coordinates. 
If Cg and have dlmensicna njj and ny^, then -I’g 
and'l'yy are 3 x ng and 3 x n;; matticcs, respective- 
ly- 


lining Eqs. 


j 


(28) , we can write 
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E ‘’“'e 


•T • 

-A -A 
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5e ^’e -e 


S “a 5a 


(w) : 


3 



uliarc 

"I ^ *'e - h " ( ^ \ ‘’'"a 

”V- '"a 

•nre nymmctrlc vomitive definite nintrlcoo of order 
itj; Olid iiy\, rei5j>ectlvely. Moreover, the nmtrlccs 
iJ,' Jjj, Ja. rE( 3 , and tag now depend on 5g and 5 a 
jlnetenU of uj; ami ua> Similarly, we liave (jjj 

il'R<Gi;). L’A " I'l: ‘ 'Ja “ ‘Ja^Sa. 

|tA)» no tliac the kinetic ciiorj;y, Eq • (9), can bo re- 
jp,nrdcd as bclnf; entirely free of spatlnl dependence. 

I 

Assuming linear elasticity, the spacecraft po- 
tuntlhl oneroy can be written in the discretized 
form 

where and Ka nic symmetric positive definite ; 
istiffncBB TnatrlcoB of order ng and nA, respective- 
|ly. Tlio roatrlccs Kj? and Ka represent intop, rals 
(Over the domain U[! and % of functions InvoivlnB 
jcpntlal derivatives of '!■£ and I'af rcspectivoly . 

Tlie iionconservatlve virtual work, Eq. (12), can 
lie discretized in n similar fashion. We shall not 
proceed with the discretization process at this 
time, but defer the question for a later soctlon. 

'i 

3. Lanranne's Equations of Motion 


Let us consider a dlBcrote (or discretized) eya- 
tcffl and denote by g(t) the conf icurntion vector of 
the entire system, 'fhen, the system l.agranBian 
can bo written in the fienerol functional form ' 


L “ T - V " L(q,q) 


(17) ; 


where it was assumed that I. docs not depend on 
time oKpllcltiy. The system admits equilibria at 
constant solutions of the equations 


3L/3g - 0 (18) 

whore 3L/3g denotes ayrabollcally a vector with the ' 
icomponents (i “ l,2,.,,,n). Without loss 

|of generality, we can n.ssume that the trivial so- 
jlutlon g ■=> 0 is n solution of Eq. (18). This la 
]i:o because one can always sliift the origin of the 
I configuration space to make it coincide with an 
equilibrium point. 

EKpandlng about the trivial solution and lin- 
earizing, the Lngranglan can be written in the 
quadratic form 

L “ q^mq + q^fq 4 ~ q'*'kq (19) 

whore n, f, and k are constant square matrices of ; 
iordcr n. Korcovor, ra and k arc symmetric. La- | 
grange’s equations of motion can bo written in the i 
symbolic form 



( 20 ) 


uberc Q is the n-d:Iraonsional generalized force vcc-' 
l.or, which can be obtained from the virtual work 
expression 

6W - g^fig (21) 


In which Is the genernUzed virtiml dluplauement 
vector. Introducing Eq, (20) Into Hq. (19), wo 
obtain the equations of motion j 

mq 4 gq 4 kq " Q (22) . 

where g - fT" - f is n skew symmetric matrix o( or- 
' dcr n. The solution of Eq^ (22) .^nn bo obtaiticd 
! in closed form, nn shown In the next section, 

6 . Modal Analygis for the llcsponiic I 


The solution of Eq, (22) can be obtained by the 
modal analysis of Itufs, 5 and 6. To this end, wc 
transform the sot of n second-order differential 
equations, Eqs. (22), into a set of 2n first-order 
equations, which , '’Mounts to working with the state 
space instead of the configuration space. 

Hence, let us define tlie 2n-dimcnsional state 
VLCtor x(t) and the associated force vector X(t) 
ns follows i 

x(t) - Ig'^(t) q'^(t)]'^ , X(t) - ig'^(t) (23) 


where 0 is the n-dlmcnslonnl null vector. Moreover, 
let us'lntroducc the 2n x 2ii matrices 



where wo note that I is syiumettlc ind C is skew 
symmetric. This permi a us to replace Eq. (22) by 


Ix(t) 4 Gx(t) “■ X(t) (25) 

I 

I Wc shall assume that Loth m and k arc positive def- 
inite, so that I Is positive definite, 

' To obtain a closed-form solution of Eq. (25), 
wc first consider the eigenvalue problem 

I 

I ' XIx 4 Gx » 9 (26) 

^ The fact that I la positive definite guarantees 
that the eigenvalues a^c pure imaginary complex 
conjugates, Ij, » im,., “ -luri ^ad that tlie as- 

sociated clgcnv^ci cu ' rre also complex conjugates, 
Xr - Vr 4 iz^, x - Iz^.. Instead of working 

with complex quaiiti i i ae , it is shown in Ref, 5 

, that the eigenvalue ,>roblcm (26) can be rcpl ced 

j by the real symmetric eigenvalue problem 

; •> ,.c «• l,2,...,n (27) 

where K ” g7^1"1g Is not only symmetric but also 

I positive definite, because I is positive definite. 

2 

! The eigenvalues Uj. of the problem (27) hove mul- 
tiplicity two. The corresponding eigenvectors are 
Y and Zf. Because I and K arc positive definite, 
the sot'of 2n aigcnvcctors and z,. (r •> 1,2,..,, 
n) arc orthogonal (with respect to'tlic matrix 1), 
They can be normalized so as to satisfy 

T T T T 

yly ‘*zlz "6 ,ylz "ziy “0 
ir is .r is rs ' ir -a >t is 

T T T T 

z Gy = -y Gz " u 6 , y Cy “ z Gz 

“b ir is -r r rs ’ is ir -s -r 

r,B - l,2,...,n 
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[ liccmilii' tl«c «i’t of vocLoni ntui In ortlio- 1 
nnrn»nl (with roiipocl. to tlio matrix 1), It conntl- < 
tiiLon n bmilti In n 2n-itli)u.>iiHlonnI vector npnr.c. 
Ilciicoi the iitntc vector x(t) con be re|U'cm>nted 
no a IJnonr combination of clieoo clgonvoctoro In 
the form 


x(t) - I + 0^(1 

r-1 


<29), 


|Whero i;i-(t) and nj.(t) ace generalized coordlnatca 
I anooclntcd with the vectoro Jfj. and Zj., renpcc- 
1 lively. Introducing Etj. (29) into Eg. (25), mul- 
1 tiplylng by mul Zp in aequunce, mid using Che 
orthogonality rDlacIoiis (28), we obtain tlic inde- 
pendent set of pairs of equations 


np(t) + 


r - 1,2, . ,n 


(30) 


where 


Yj.(t) - y^Xd) , Zj.(t) » Jx(t) 

r " 1,2, . . . ,n (31) 

arc generalized forces associated with the gener- 
alized coordinates C;[^(t) and t)p(c) respoctivcly . 

Equations (30) can be solved for the pair of 
generalized coordinates fp(t) and n_(t) indepen- 
dently of any other pair. Introducing this solu- 
tion into Eq. (29), wo obtain the complete res- 
ponse 


x(t) 


n ft „ 

I t J t(y^yj 

r"l ^ o 


+ Z^Z^)X(t)C 0B (Uj.(t-T) 


^Xr5r *■ “p(*^"f)]dT 

+ (y_y^ + z_2^)Ix(0)cos m t 

•’TtV c 

+ (Xr?p " Zj.’y^)lx(0)sln tU|.t} 

where x(0) Jyj the initial state vector. 


(32) 


The decoupling procedure can bo written in ma- 
trix form. To this end, let us introduce the 2n- | 
dimensional generalized coordinate vector : 


^n V 


w - ni 1,2 H 2 
as well as the corresponding modal motrlx 
P 


f?l 5l ?2 52 •' 


y z ] 
in 'n^ 


BO that Eq. (29) can be written In Che compact 
form 


Pw 


(35) 


Horoovor, the orthogonality relations con be com- 
bined into 


T 

P IP 


*2n 


(36) 


wbcrc l 2 p is the unit matrix of order 2n. Pre- 
tmiltlplylng both sides of Eq. (35) by P’'^I ond 
using Eq. (36), we conclude chat 


w - l'’’lx 


(37) 


Next, let us define o block diagonal matrix 
A as foilowa; 


1 

0 

— 

0 ~ 

0 

“2 ^ 

— 

0 










0 

0 

— - 

“n \ 


» i " 


0 

-1 


(38) 


With tills notation, Eqs. (30) and (31) can be 
written In tlio compact f jrm 

I w " Aw + » X (39) 

7. Rosponso of Ui on trolled System to 
Imnulsiyc Excitation 

bet uB consider the case in which system (25) 
is subjected to on impulslvo force at time 0+ 
while in equilibrium, x(0) “ 0. The force vec- 
tor can be written in the form 


X(t) •> X6(t) 


(AO) 


! where X is tlic magnitude of the impulsive force. 
Introducing Eq. (AO) into solution (A9), wo ob- 
tain 

<(t) - I [ + 5r5r)5A(T)sos Uj.(t-T) 

r"l 

+ (y^z^ - Z|.y^)X6(T)sin iiJj.(t-T) )di 


I [(y Xp + Z Zp)Xcos (ii t 
r»l 

^^r^r “ 5r?r^ 


(Al) 


wliich sliDws tliat the motion of the system consists 
of a superposition of harmonic motions at the nat- 
ural frequencies Up. 

Tlic Impulsive force X6(t) can bo shown to 
cause n motion analogous to that caused by an Ini- 
tial excitation. Indeed, let us introduce Eq. 

(AO) into Eq. (Al) and write 


(33), 

I 

(3A)' 


lx(t) + Gs(t) - X6(t) 


(A2) 


betting tile durotlcm of the impulse be At end 
integrating Eq. (A2) with respect to time, we ob- 
tain 


fAt . 

j X(t)dt + G 


fAt « fAt 

x(t)dt « X A(t)dt 

» o * n 


(A3) 


For small At, the second Integral on the left 
side of Eq. (A3) Is negligible, so that 


fAt , 

ijn I x(t)dt - l|x(0+) - x(0)) - Ix(Ol) “ X 
•*0 ' o 

(AA) 

From wlilch we conclude that the Impulsive force 


lijn 
At- 
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|iro(Uicuii Llin ci|ulviilvnt Initinl oxcitntlon 

x(OH) - r^X (A5) 

It cnii be ennily vertCloil tiint 1C Cbo Initial ex- 
citation (45) in Innurtud Into K<]i (32) Instead 
oC the Corce (/«0) thn result would be tlic samci ; 

I 

8. System with Proportional Control | 

Lot us ooooino that the iiyntoin under considera- 
tion was subjected to an impulsive Corce rcsiiltlnn 
In oscillation according to Sqi (^il) . Let us ! 
Curtber ossunic that the response exceeds a given 
amplitude. Itecause the oscillation persists with- ' 
out attenuation) the response must bo considered { 
unuatlsCactory I so that It Is deemed necessary to 
attenuate It by means o£ active controls. In this 
section) we stinll use proportional controls and 
In the next wo shall use on-oCC controls. ' | 

Denoting by U the control vector) the system 
dlCCcrcntlal equations oC motion can be written ! 
in the Conn 


Combining Hqu> (AO), (A9)i and (S2) , we can write 


Yp(t) - -cCj(t) - -c^^lx(t) - y’^U(t) 

1 

1 

z (t) •< -cn_(t) » -czjxx(t) " z][u(t) 

I 

i 

1 

L' 1 , 2 , . . , , n 

(53) 

from which we conclude that 

1 

i 

U(t) " -clx(t) 

(5A) ■ 


! 

or, the control vector U(t) is proportional to 
the vector lx(t) . ' 1 


9 , System with On-OCC Control 

Proportional control lias one drowbnek) namely) 
it must opernte continuously. A scheme without 
I this drawback lo on-oC£ control. The control law 
assumes n region o£ deadband bared on the recog- 
nition that wltliln some tolerance small oscllln- 
> tlons arc acceptable. 


lx(t) + Gx(t) » y(t) (A6) 

which is subject to the initial conditions x(0) . 
Using the approach oC See. 0) we can reduce' the 
simultaneous set (A6) to the independent set 

Cj(t) - w^y^d) - Yj.(t) 

r « 1,2) ... ,n (A7) 

yj.(t) + Wj.5j.(t) “ Zj.(t) 

wlierc 


I Lee us consider once again the system (A6) 
i and the decoupling procedure (A7) and (A8), but, 
' by contrast, we assume a control In the form 


U 


n 

I 

0“1 


“s “ 


(55) 


where Ug Is a nonlinear function of gg to bf. 
specified sliortly. Introducing 8q. (55) Into 
Eqs. (A8) and considering the orthonormallty re- 
lations (28) ) wo conclude Chat 


Yj.(t) “ )!^l|(t) ) Z^(t) “ z^y(t), r " l,2,,..,n 

(A8) 

Next, let us assume proportional controls 
Yj.(t) ° -cK^(.t) , Zj.(t) ° -cn^(t) , 

r »• 1,2, . . . ,n (A9) 

so that Eqs. (/,7) can be rewritten in the form 

Cp(t) - w^n^(t) + " 0 

r=l,2,...,n (50) 

nj.(t) + u^Cj.(t) + cnj.(t) = 0 

The solution of Eqs. (50), obtained by the La- 
place transform method, Is 

5^(t) “ c”‘^‘^U;^(0)co3 u„t + nj.(0)sln Wj,t] 

i'ij,(L) •* e”'^*^(-f;^(0)sln u^c + nj.(0)cos u^t] 

i 

j r-l,2,,..,n (51) 

I no that the response dies out with time. 


- 1 T 

Y^<t) •‘I — y i 2 .u - 0 
r w ir -a s 


G"1 B 

n 


Z(t)" y — kIzu “ — u 
r' ' w <r >0 8 u r 

o“l B r 

so that Eqs. (A7) reduce to 


(56) 


- 0 


, t r “ (57) 

i ~ “r " ° 

i 

I Next, let us specify that Che function u,. la 
1 given explicitly by 


Uj. “ ( 0 , -d^ <. Ij. <. 

( K . 


(58) 


n < -d 
'r r 


lienee, the solution of Eqs. (57) must be obtained 
acpnratcly for the three intervals Indicated above: 

1. For |n^! S. ‘'f* •‘‘I-'’* (57) reduce to ! 


Thu question ruraaJiis ns to how the proportional . 
control on the decoupled coordinates l|-(t) and | 
i|,.(t) rcl.ntc to that of the st.sue vector x(t). To | 
1 thin end, we multiply Eq. (29) by and'zel, ir. 

I sequence, coiusider Eqs. (28) and obtain 

lc-(i-) " y3lx(c), II (t) «• z^Ix(t), r “ 

^ ___ I ( 52) 


i^(t) - Uj.nj.(t) - 0 

r - 1,2 n (59) 

nj.(t) + " ° 

which can be sliown to have the solution 
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c^(t) - Cp<0)con Wj.c + n^<0)nln 

ilp(t) “■ -Cj.(0)oln Uj.t + nj.(0)coit 

r-l,2,...,n (60) 

il. For > dp. Up « -kp, Eqn. (57) bccono 

Cp(t) - Upiip(t) - 0 

, k_ r - 1,2,., ,n (61) 

tip(t) + UpCp(t) + ~ " 0 

whlcli have die oolutlon 


VO- ^ + 
r 


(^(0) 


_r 

2 

“rj 


cos u t + n (0)6in u C 
r r r 1 


tlj(t) 


Cp(0)+-| 

"r_J 


sin w„t + n„(0)coo w t 
jf r r 


r « l,2,.,.,n (62) 

ill. For Hj. < “dp. Up •* kp the response is ob- 
tained by simply rcploclno "kr by Hkr In 
Eqs. (62). lienee, 


Cp(0^^ 


c,(0)--| 

“r. 


COB u C + n (O)sln h) t 
r r r 


<1^(0 


- 


_r 

■L 


sin Dpt + tip(0)coa Upt 


1 , 2 ,.. 


(63) 


l"i 


The behavior o£ the solution, Eqs. (60), (62), 
and (63), can bo discussed most conveniently in 
the phase plane n,- vs. f.r. From Eqs. (60) , we { 

' ni 7 

I conclude that If /(^(O) + n,-(0) < d^, then the 

I trajectories represent circles with the centers ^ 

|nt the origin and with rodll + rif(O). I£ 

idle motion Is Initiated In the region Pj. > dp, then 
tlrora Kqs. (62) we conclude that tlie trajectories 
o re clrclcg centered ot S r " “kr/Uj end with radii 

[f,r(0) kf/uj ] + n,;(0). On the otlier hand, 

(rom Eqs. (63) we conclude that the trajectories 
Initiated in the region Pp < dp are circles cen- 
tered at f,p «■ kp/u“ and with radii . 

/[(r (0) kp/u)‘^]2 + n^(0) . A given motion Inl- 

! tinted at some point Pp(0) > dp, (p(0) < 0 will 
liollow a circular trajectory until It reaches the 
I liorlzontal line Pp “ ilp, when the control Is re- 
{moved. If at this point r.p < 0, then tlie trnjec- 
j lory will tend to move clockwise on a circle cen- 
tered at the origin, which will take the motion i 

hack Into the region Pp > dp causing the control ' 

I to he actuated ng.iln, Kepctltlon of this motion i 
I pattern results in chattering along the line Hp “ I 
Idp. If Che trajectory Initiated at nr(0) > dp j 
jhlts the line Op “ dp at a point for which (p > 0, 
j then the motion will continue on a circle with Che ' 
center at the origin until It reaches Pp = -dp, so 
that now chattering occurs along the line Pp ° 
l-dj.. Figure 2 shows these trajectories along with 
I^DOmc other iiosslhle c.sscs. 


To prevent chntcnrlnit, one may wish to delay 
the time no that tlio cuiitrolii aru removed ulien 
tho trajcctorlcn arc Iniildo the deadband Interval, 
This in oqulvnlcnt to Introdticlnn a phaiia iingli) 
in the iioluUona (62) and (63). Figure 2 nliown in 
dashed line tlie syoccni behavior Cor a time delay 
corresponding to n pliaoo lag of 10° , 


Tlie on-off control wltli dendbund and time de- 
lay corrcii]iondlng to n decoupled mode la llluntrated 
in the block diagram ahown In Fig. 3. 

It nliould be pointed out Chat the abovo anniy- 
olo, including tlie phase plane representation of 
: the motion, would not have been possible without 
I tho decoupling procedure, ! 

10 . Reconstruction of tho State Vector 
from Available Outnucs 

1 1 

Regarding the control as nn external excitation 
and using die analogy with Eq. (61), Che equations 
. of motion of a controlled spacecraft subjected to 
1 external excitation can bo written in the form 


lx + Gx “ X + U 


(66) 


Unlike the force vector X, however, which depends 
on time alone, the control vector U depends spe- 
cifically on the state variables, 'in Sec. Q, we 
otudlcd two cauQU, namely, tiint in v.Ulch U is n 
linear function of the state Vector and that In 
which y is n nonilnoar function, liccauso now tho 
Input to tlie system is X + U instead of X, the 
decoupled equations of motion can bo written in 
tho symbolic form 


w - Aw + 1>‘ (X + y) 


(65) 


The above control Is predicated upon the knowl- 
edge of tho stote vector x(t). Quite often, how- 
ever, the state vector Is'not completely known, 

BO that 0 method for its estimation is highly de- 
sirable. Such a method uses anotiicr dynamical 
system known as an obBorver . 

The dlscusBion of tlie observer can be conven- 
iently presented In terms of tho uncoupled system. 
Lot Wp, denote a vector of monsuremonts correspond- 
ing to the time derivative of the uncoupled utato 
vector w. The object la to canntruct nn observer 
capable of yielding a good estimate of the state 
vector. Such nn observer should be a dynamical 
system resembling the dynamical system (65) and 
should depend both on the Input X + U and the 
measurement w^,. In general the measurements need 
not be the complete set of states. The observer 
approach can easily be extended to include the 
effect of measurement errors thereby resulting in 
a filter approach. Ilcncc, let us consider nn ob- 
server described by the following vector differen- 
tial equation 


w - A„w + II Cl + N„P‘(X + U) 

- Q' O'Cl o - - 


(06) 


I where w denotes the observer-constructed state 
I vcctor'nnd Aq, 1!q, and Nq ore block-diagonal mn- 
I trices to bo dotcrralued so that the observer ex- 
I mbits the desired behavior. Moreover, the mea- 
surements vcctur w„ Is related to the stare voc- 



tor w by 



w - C„w 
-ra 0* 

(67) 
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iwlii'io To III iu'iii''iil ly I* I'lH'lmii’.uliir mtii rlx ulitcli 1 

lU'ID'iiilii uii Llif I lhii iii;li'rliiLliu 1)1 Lliii tiii'iiiiiirlnn 

iluviccti. intrixluL'lii}’ IlitUi (G!i) mill (G7) liiLo 
l'i|. <C6) iiud {)ubLr<)ctl.tii’, tlia ruuuit from Bq, (65) > 
uc obtnttt 

w - w - (i - !loCy)Au - AqM + a - CqHq 

-NQ)r^(X + U) (60) 

wlipro 1 ill tlio IdentlLy matrix, boctini! thci mn- 
trJcoB Aq, Uq) nml Mq imtiofy ttic uqiiiif ivmB 


0 - Cq»o)A h Aq 


^ ■’ °o"o “ "o 


(69) 


Eq> (68) rciiucco to 

w - w •• Aq(w - w) (70) 

Tho obiicrvcv (lynnmics is slmuiatctl in Flgi A< 

Next, let us Introduce the notation 

c(t) “ w(£) - w (71) 

where c(t) is known ns tlio cstlnintion error . i,c., 
the diCCcrunce between the actual state vector and 
the reconstructed state vector. Introduclne Eq. 
(71) Into Eq. (70) I we obtain 


c(t) " AQt(t) 


(72), 


no that if the elijenvnluoa of Ag have negative reol 
parts, tho error decays witli time. There are no 
aiipnrcnt restrictions on tho eigenvalues of Ag un- 
til iticasurcmcnt errors ore taken into account. 
Typically, the cigenvaluca of Ag would bo "Coster" 
than the syntciii eigenvalues by a factor of 5 to 10. 

I 

I 

Ccnerolly, one chooses the nntrlx ilg so thot 
tho matrix Ag has the desired eigenvalues. In par- 
ticular, Uq is ctioncn ns n block diagonal matrix. 

We recall that the matrix A is itself block-diago- 
nal, os can be seen from Eq. (38). j 

Next, let us nsnumc that some components of tho 
state vector x(t) have been measured by means of 
on-bonrd sensors, such ns rate gyros, occolerom- 
cters, etc. Denoting the measured state vector 
rate by ^(f)) where 


-m 




we can use Hq. (37) and w'ite 




l-^lx^(t) 


(73) 


(7A) 


lAssumlng for simplicity chat Cq is tt)c identity 
'wnlrlx, which impJles complota observability, and 
introducing Eqs, (69) and (7A) into Eq. (66), we 
oh till 11 the. observer equation 

• 

w " A_w + B„l>'*^Ix + (1 - B,.)l'^(X + U) 

0- 0 -Itl Q - • 


Aq id the dXnf^onnif i*c. > 

Aq " dlag(aj^ Oj, 0J. .. 

1 tbVlowE immediately that 
-1 


1 - K 


A^A 


“ blork-dlag 


(X + y) 

(75) 

so that the 

matrix 


(76) 

n n 

0 “U,./'*'.. 



(77n) 


lljj » block-ding 


r ^ 


Introducing tlio notation 




whore 


9 " ••• 

nnidoring Eq 
form 

A * 

C " o C + Q, 
'r r^r ^(r 


‘>f.n 


(77b) 


(78) 


(79) 


and connidoring Eq. (76), Eq. (75) enn bo written 
in tho form 


r " l,2,i«*,n 


(80) 


\ " ^'’r 


: which lino tlie general solution 

I M * _ ft 


C^Ct) 


tp(0)o' 




Hp(t> ■ + 


'0 

f ont(t-T)* 

* Q 


Q^,(T)dT 


1|2.. 


(81) 


Choosing Of nnd 8^ (r ° l,2,.,.,n) an complex 
numbers with negative reel parts, we conclude from 
Eq. (72) that tlio error c(t) reduces to aero with 
time, so that the decoupled obuorver state vector 
w can be used to determine the behavior of tlio de- 
coupled system state vector Note that to ob- 
tain tho actual observer stote vector we can write 


x(t) - I’w(t) 


(82) 


The vector Q consists of one part due to the 
measured state'vcctor and external disturbances 
another part due to tlie control 


where 


Q + Q 
3x 3U 




(1 - Bq)P*U 


Bq)p'‘^X 


(83) 


(8/i) 


Hence, the classical separation nrlnciplo Is illus- 
trated In that the obscrver^can be used to actual- 
I ly control tho system. If gi, is taken to slrau- 
i late the proportional control of See. 8 or the on- 
off control of See. 9, then the actual control 
vector is obtained from the second of Eqs. (8A) in 
the form 


U 


ip(i - BJ-ig., 


where 


(1 




-1 


- block-ding 


0 

_-m /a 0 _ 


(85) 


( 86 ) 


11. Control of a Spinninp, Flexible Siincecrntt 

Before proceeding with the derivation of Ln- 
IgrnpGc'.n equations of motion it will prove con--- - 
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EI^’HODllClBiU ' 

nRK.'fiSl.AL FAVii' 


lo iiti'tUlfy thi> Hyiili'ni ni'iti’riil l/.i'il coor- ] 

I AioMinil IliiiL llu* lull 1 1 1 li* tmivi'it In il 

jcirculnr orhlt nround tlio cnrtli iiml timt tlio cun- 
(lor of ntiMio of Llu> miU'lUU> colnclilon wUli tUo 
ici>ni.(’r of riiiiiti I' of tlio plalfona fov nil prncticnl 
I pill poitoii , ihi' vc'clor V|)p ciiii III' coon 111 vi'L'il no boloH 
known. Moroovor, wo oli.iU nnminu> tlmt lUoro nro 
inn lac’niliiTii rotnllnft relntlvi* to tlin plntforn mill 
I Hint tlioro nro tvo I'JniiLlc piinolH nitnclicil to thn 
j'lntform, Tlinim pnm>ln nro iiymmctrlc with rciipcct 
I til llii' point I' omi thi'lr niotlun In notliiyimiu'trlc. 

Ill follown tImt nil qunnticii'u with inibncriptn It 
tnml A cmi bo jRiiori'il in thu klnotlc oiusrp,y Kq. (9). 

I in ndiUtioiii wc hnvu 

! "’/S'K " ^EC " ° ' h " ° 

CoiinJdorinii thu nbovc noBiimptionni no well no Eqo. 
(lA), and Icnorliifi conotont tormni Kqa, (9) roducun 
to 

T " 2 ?p 5p ^ 5 e ^’e 5e ** 2 p **PE *B 

I Next, let uc aBBUOG that In equilibrium tUo 
iplntform axca Xpyj,Kp rotnto rein live to the inor- 
itial epnee XYZ with tho uniform nneuinr velocity 
Ifl nbout 2 j>, where lo pnrollel to Z. Then, if 
luc cotiBlder n set of atixllinry nxao xptypizp' ro- 
ItntliiB rolotivo to the Inortlnl npaco with tho nn- 
lEiiJar velocity si nbout Zni. whore Zpi Is pnrnllcl 
I to Z, and if we denote the nn'n’.'i ir velocity of 
Inxea XpypZp relative to xpiyp y., i by wpi tlion the 
i nneuinr velocity of the platform axes Xp/ptp rela- 
tive to tho inertial apace can be written in the 
form j 

5j, “ nj + Wp 

wbere 1 la the vector of direction cosinca between 
axia Zpi nnd axca XpypZp. AssuralnB that nxea Xpyp 
Zp nro ohLalncd from nxea XpiypiZpi by monna of 
tlicn rotations U2 about yp, 0 ^ about Xp, nnd O3 
jnbotit Xp, in that order, then j 

I >p 

» [nOj^nO^-cO^BO^cOj aO j^cOj+cO^aO^sOj cOj^cOj] 

(90) 


nnd 



CO3 SO3BO3 0 ~ 




i’p " 

-BO3 CO3CO3 0 


°2 

( 91 ) 


0 -aOj^ 1_ 


-° 3 - 



whore aOj - nin 0 ^^, cO^ “ cos 0 ^ (1 *• 1 , 2 , 3), 

heS'reiieutinB the clnatlc diapluccinontu of the 
pnncls as folloua: 


“E “ ° • ''e " ^l*^! ’ ’'e “ ‘^a '^2 '^3‘^3 

Iwhere ijij la ‘he first In-plniie mode, l ^2 is tho 
I firm our-of-plnne mode, and i('3 in tlio first tor- 1 
' iiionni mode about x,,, iiitroilucIiiB Kqa. ( 89 }-( 92 ) | 

Unto Eq. ( 88 ), and Ilnoarlzlnp,, wc obtain 

T •> I n^tCB - C)D^ + (A - OOj - 2a0j^i:3 + I 


+ m/.Il + O[-a 6^02 + (11 - j 

- "°2 ■’* ' 2 '^Vl 

I 

I 

+ EtjCj + iBjCj) - bO^Cj nOjCj^ + eOj^ij (93) 

where A, II, C nro tho momentu of inertia of tho 
cntlru undeformed npneccrnft nbout nxoii xp, yp, zp, 
ruupnctively, nnd . 


" 1 ’'e^ 1 ‘’"e ' " 1 *E* 2 '‘"’k • 

" I ^e*d‘''^ ’ '”t " I ^ 

• «!._ • fn. . 


(9A) 


“E 

I In termu of tho modoo Indicated by Eqa, (92), the 
potential enorcy Is 


V»i I m.Afci 

' i“l ^ ^ 


(95) 


where Ai (1 ■ 1,2,3) nro tho natural froquoncloa 
of the panel nssoclntcd with tho modes (1 ■> 1 , 
2,3). 

The nonconsorvntlvo virtual work enn be written 
in terms of the gonarnllzcd forces nnd virtual dls- 
placcmonto ns follows: 

I “ I ^'=i> 

I i«»l 

where nnd F 51 (1 - 1,2,3) nro tho Bunornlized 
forces. 


LnBrnnBc'o equations of motion can be written 
in the B°ncrol form 


d_ ( ^h\ ah 

iL ^ ah \ ai. 

VcJ “ “^1 


« F 


01 


1 ” 1,2,3 


(97) 


'cl 


From the third of Eqs. (97), we obaerve that If 
F „3 “ 0, then O 3 is an iBnorabie coordinate, ao 
tliut 


ai 


— " CO 3 + ■> (I " const 


, ao 

! J 

I 

, lienee, introduciit); 

8 - aq^ 


( 98 ) 


( 99 ) 


j into tlic kinetic unerp.y, wo obtain j 

T - “ n^Ub - C)0^ + (A - OO 2 - 2e0j^53 + 2b02q2 

i 

+ Pij^cJl + nt-AO^O^ + (8 - 082 °! " *’^°l '‘2 ® 1 ^ 2 ^ 

2 ^ 

" ®° 2 '^ 3 ^ 2^*°1 f Knij^ - ; 


t 
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+ ■' "'353I ■ '‘**2^2 


( 100 ) 


Sicxtt let lilt Intraduca ttio configurntion vector 


3 - lOj Oj r, C3I 


( 101 ) 


mill ihu iimiocinU'd mmcotincivdtiVQ force vector 

T 


9 “ ^^'ol '02 f’u ^2 •'cdI 


( 102 ) 


Then I.nr,rmt|',c'H uquntlonHt ilq. (20), nooumett tlio 
foru (22) D In wliicli 


0 

II 

0 

-l> 

0 


0 

-b 


D 


1 2C 
0 e>2 

0 0 


'3J 


0 >06/1 7 red 11 


A.J 


2 


A3 >• 0.0227 red ti 


0,jl7A3 


red (I 


Siniilntlonn were mmlo of the nyiitein reeponne fort 
[D tincontrollod iipncocrnft, 2) p report I oiml control, 
1 3) on-off control with ilQiidbiind. Fiijurcn An, 4I>, 
pn, Sli, imd On, fib nhow typlcnl compntor plotn for 
the nntntJon nnttio Oj imil tliu unt-of-plniu> mode 
I for tlio tliroo cmiou, rciipuctlvoly . Tho Inttlnl 
condltlono nnd tlio vnrloiiii control pnramotorii woroi 

! 0^(0) - 02(0) - 10“''* rnd, Cj^fO) - ^^(O) - 1:3(0) 

I » lO"^* n , c •• 0.1 n”^ 

d, “ d, - 10"^ , d, - d, - d, - 10"'' ' 

^ ^ -3 3-2^' 5 -3 2-2 

I k, •• 2 X 10 ■’ uf 0 , k, - 2 X 10 u, B 

-A 2 -i ^ -A 2 -2 

, “3 “ 10 >3 " : • ^ - 10 “A “ • 

I kj •• 10" ' u‘ , 

11).' Incliinlon of Plp.n. A, S, fi In me roly to iiliow 
typical rcBultii. I'lio rcnultu nrc not meant to ru- 
prenent optimal control. Indeed, Llie on-off con- 
trol rcmilta could bo r.rcntly Improved by a reduc- 
tion of the dead band constants. 


0 C-A-U 000 


12. Concluolonn 


r. » R 


(C-A-B) 

0 

0 

0 


0 0 0 -e 

0 0 0 0 

0 0 0 0 

c 0 0 0_ 


(103) 


C-ll 0 


0 0 o 


0 C-A 


0 -b 0 


k " n 


0 0 

0 -b 
o 0 


la^(Aj^'-l) 0 0 

0 "'2^2^ ® 

0 0 mjA* 


I where A* ■ Ai/tl (i “ 1,2,3). 

Ui! obtterve that the equation for the coordinate 
in Independent of the other oncii, bo that Its 
I noltitlon c.nn be ebcnlncd Independently. For olin- 
tpllclty of computer proprnmmitir, , however, we choose 
'not to treat iicparntely but ns port of tho for- 
□ul.ntlon (103). I 

i t 

I The Boliitloii of Eq. (22) with nnd without con- : 
I trolo and with m, i,, and k .in Riven by Eqa. (103) 

' lollcrws tiic p.nttern eBtnbllcihcd In Sees. 6-10. 

I 

The above formulation was used to detcimlne the . 
renponne of n spacecraft with the following parnra- 
ctoras 

A - 1,000 kg m^, B = 6,000 kg 
C •» 8,000 kg m^, 0 “ 0.6 rnd b“1 

Tho pnncln wore modelled by the finite element 
'method. The first natural frcq"oncicn for In- 
|pl.3nc, out-oC-plnne, and tornlonal vibration are 


Tills paper develops o modal procedure Cor the 
control of n flexible spacecraft exhibiting gyro- 
scopic behavior. Control via decoupling has dis- 
tinct computational advantages over control of the 
:couplod oyotem, particularly frr large order syo- 
I terns, as It pcnnlts the use of methods of solution 
jgcnernlly associated with second-order systemo. 
Design procoduroii nro domonotrntod for two types 
|oC control slgorlthma, propurtlunul control and on- 
;of£ control with dead band, 

I 
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